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Research  has  been  carried  out  in  two  areas:  first,  the  extension  of 
Pisarenko’s  method  to  multidimensional  (m-D)  signals  and  arrays;  second, 
the  extension  of  the  maximum  entropy  method  (MEM)  to  non-uniformly 
sampled  ra-D  signals  and  arrays.  Most  of  this  work  has  been  of  a  theoretical 
nature. 

A  procedure  for  calculating  Pisarenko’s  estimate  has  been  derived.  The 
form  of  the  estimate  is  that  of  sinusoids  (or  plane  waves)  in  a  background  of 
noise  with  known  spectral  content,  but  unknown  total  power.  The  problem 
can  be  reduced  to  a  linear  program  that  is  computationally  similar  to  the  one 
used  for  designing  two-dimensional  finite  impulse  response  (FIR)  digital 
filters. 

In  the  process  of  deriving  Pisarenko’s  spectral  estimate,  one  is  led  natur¬ 
ally  to  consider  the  "extendibility"  problem:  Given  a  finite  set  of  samples  of  a 
presumed  autocorrelation  function  (acf),  does  there  exist  a  positive  power 
spectrum  that  exactly  matches  the  given  acf  samples?  This  extendibility  prob¬ 
lem  is  important  in  practice.  Often  the  autocorrelation  function  is  estimated 
from  data  samples  and  it  is  necessary  to  have  a  test  to  determine  whether 
such  acf  estimates  are  extendible.  One  possible  test  is  to  compute  Pisarenko’s 
estimate  under  the  assumption  that  the  background  noise  is  white.  In  the 
course  of  the  computation,  the  white  noise  power  will  be  determined;  if  it 
turns  out  to  be  negative,  the  acf  estimates  are  not  extendible,  if  positive  they 
are  extendible.  If  the  noise  power  were  zero,  then  the  power  spectrum  would 
consist  only  of  impulses.  Improvements  on  this  test  will  be  needed  in  the 
future. 

The  papers  by  Lang  and  McClellan  [1,2,3]  develop  the  theory  of 
Pisarenko’s  method  in  detail.  Connections  to  some  practical  array  processing 
problems  are  also  given.  The  review  papers  [4,5]  by  McClellan  also  discuss 
Pisarenko’s  method  in  the  context  of  the  general  multidimensional  spectral 
estimation  problem. 

The  second  major  topic  of  research  was  the  development  of  a  general 
multidimensional  MEM  spectral  estimation  algorithm.  A  number  of  such 
algorithms  have  been  proposed,  but  the  research  under  this  contract  [6,7,8, 9] 
is  unique  in  that  a  rigorous  mathematical  approach  was  used  to  exploit  fully 
the  convex  nature  of  the  MEM  optimization  problem.  In  particular,  a  general 
necessary  and  sufficient  condition  for  the  existence  of  the  m-D  MEM  spectral 
estimate  has  been  obtained.  This  existence  condition  comes  about  when  a 
dual  entropy  minimization  problem  is  derived.  The  dual  problem  is  particu¬ 
larly  important  because  it  is  a  finite-dimensional  convex  optimization  prob¬ 
lem;  whereas,  the  original  MEM  problem  is  infinite-dimensional.  Much  of 
the  previous  work  on  MEM  used  some  form  of  the  dual  problem,  but  not  the 
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form  derived  in  this  work.  A  most  important  point  is  the  convex  nature  of 
the  problem  which  guarantees  that  traditional  numerical  optimization 
methods  such  as  steepest  descent  or  quasi-Newton  will  converge. 

A  final  area  of  investigation  is  the  use  of  the  new  m-D  spectral  estima¬ 
tion  techniques  in  a  practical  setting.  Some  work  [10,11]  has  been  published 
on  the  comparison  of  MEM  with  interpolation  for  bearing  estimation  with 
non-uniform  arrays.  This  research  has  shown  the  superiority  of  the  general 
MEM  algorithm  for  synthetic  signals,  consisting  of  sinusoids  in  noise.  Work 
is  presently  underway  to  apply  these  algorithms  to  data  from  seismic  arrays. 
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